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Abstract 

We  consider  the  reliability  of  a  single-unit  system  whose  cumulative  damage  over  time  is  a  continuous  wear  process 
{X(t):  t  0}  that  depends  on  an  external  environment  process  {Z(t)\  t  ^  0}.  We  explicitly  derive  the  failure  time  distribution 
and  moments  in  terms  of  Laplace-Stieltjes  transforms  by  analyzing  the  Markov  additive  process  {(X(t),Z(t)):  t  ^  0}  and 
demonstrate  its  applicability  on  an  example  problem. 
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1.  Introduction 

In  this  paper,  we  investigate  a  single-unit  system 
that  accumulates  wear  over  time  due  to  the  influence 
of  a  random  environment  that  may  be  modeled  as  a 
continuous  time  stochastic  process  {Z(t):  t  ^  0}.  In 
particular,  we  assume  that  the  wear  rate  of  the  system 
depends  explicitly  on  the  state  of  its  random  environ¬ 
ment.  Assuming  that  repairs  do  not  take  place  to  re¬ 
store  the  condition  of  the  system,  the  cumulative  wear 
up  to  time  t  may  be  characterized  by  a  nondecreas¬ 
ing  stochastic  process  {X(t):  t  ^  0}.  We  discuss  the 
properties  of  this  process  in  Section  2. 

The  system  begins  its  lifetime  in  perfect  working 
order  but  experiences  wear  under  the  influence  of  its 
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random  environment  until  the  state  of  the  system 
exceeds  a  fixed  threshold  value  x,  at  which  time  it 
fails.  We  denote  the  time  until  failure  by  a  random 
variable  Tx  and  derive  explicit  transform  expressions 
for  the  cumulative  distribution  function  and  the  nth 
moment  (n  ^  1 )  of  this  random  time  when  the  envi¬ 
ronment  process  is  assumed  to  be  a  temporally  homo¬ 
geneous,  finite-state  Markov  process.  Our  main  results 
are  obtained  by  demonstrating  that  the  joint  probabil¬ 
ity  distribution  of  the  Markov  additive  process  (X,Z), 
conditioned  upon  the  initial  state  of  the  environment, 
satisfies  a  first-order  linear  partial  differential  equa¬ 
tion  (PDE)  which  is  directly  solved  via  transform 
methods.  Next,  using  the  derived  distribution  of  the 
failure  time,  we  characterize  the  Laplace-Stieltjes 
transform  of  the  failure  time  moments.  The  results  of 
this  paper  are  well-suited  for  practical  implementation 
with  existing  single-  and  multi-dimensional  Laplace 
transform  inversion  algorithms. 

The  modeling  of  stochastically  deteriorating, 
single-unit  systems  has  been  studied  extensively  in 
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the  literature  by  a  number  of  researchers.  The  primary 
approach  has  been  that  of  stochastic  shock  and  wear 
models.  Shock  models  can  generally  be  described 
as  follows.  A  system  is  subject  to  shocks  of  random 
magnitude  at  random  time  intervals.  The  cumulative 
wear  to  the  system  is  characterized  by  the  sum  of 
discrete  shock  magnitudes  accumulated  up  to  some 
time  t.  These  models  have  been  primarily  used  for 
prescribing  preventative  maintenance  procedures,  e.g. 
optimal  replacement  policies.  For  an  extensive  sum¬ 
mary  of  papers  addressing  this  issue,  the  reader  is 
referred  to  the  excellent  survey  by  Valdez-Flores  and 
Feldman  [16].  Gottlieb  [6]  gave  sufficient  conditions 
for  the  shock  process  in  order  for  the  system  failure 
time  distribution  to  have  an  increasing  failure  rate 
(IFR).  Waldmann  [17]  considered  a  shock  model  re¬ 
lated  to  an  environment  process  that  describes  varying 
external  factors  and  found  an  optimal  replacement 
policy  under  this  assumption.  More  general  Markov 
renewal  shock  models  have  been  analyzed  apart  from 
the  optimal  replacement  problem.  For  example,  the 
probability  distribution  of  failure  time  for  two  distinct 
Markov  renewal  models  were  considered  by  Shan- 
thikumar  and  Sumita  [13]  and  later  extended  by  Igaki 
et  al.  [7].  Some  examples  of  shock  models  influenc¬ 
ing  multi-unit  systems  include  the  papers  by  Ride 
[11],  Nakagawa  [10],  and  more  recently,  Stadje  and 
Zuckerman  [15], 

The  seminar  paper  by  Esary  et  al.  [5]  considered 
processes  subject  to  continuous  wear,  relaxing  the 
shock  model  assumption  that  deterioration  occurs  only 
at  discrete  time  points  and  in  discrete  magnitudes. 
When  the  continuous  wear  process  X  is  assumed  to 
be  a  nondecreasing  Markov  process  whose  resistance 
to  wear  decreases  with  time  and  wear,  it  was  shown 
that  Tx  has  an  increasing  hazard  rate  average  (1HRA) 
distribution,  f'inlar  [4]  generalized  most  of  the  mod¬ 
els  of  [5]  by  demonstrating  that  the  process  (X,Z)  is 
a  Markov  additive  process  and  gave  several  examples 
of  such.  The  first  example  relevant  to  this  work  (cf. 
[4],  p.  199),  considered  the  case  when  Z  is  a  general 
Markov  process  with  a  finite-state  space  and  X  was 
assumed  to  increase  as  a  Levy  process.  Addition¬ 
ally,  random  shocks  were  assumed  to  occur  at  each 
environment  transition.  In  such  case,  the  X  process, 
conditioned  on  Z  follows  the  same  probability  law  as 
an  increasing  Levy  process.  The  second  example  (cf. 
[4],  pp.  201-202)  is  similar  to  the  problem  discussed 


herein  where  A  is  a  continuous  additive  functional  of 
Z  and  the  failure  time  is  a  first  passage  time  for  the 
X  process.  More  recently,  Abdel-Hameed  [3]  exam¬ 
ined  properties  of  the  time  to  failure  when  the  wear 
process  X  was  assumed  to  be  an  increasing  Levy 
process.  Singpurwalla  [14]  gave  an  excellent  sum¬ 
mary  of  a  variety  of  stochastic  failure  time  models 
for  systems  in  a  random  environment  and  particularly 
noted  the  difficulty  of  implementing  many  of  these  in 
a  practical  setting.  Liminos  and  Oprisan  [8]  provided 
mathematically  rigorous  results  for  several  reliability 
measures  for  semi-Markov  systems;  however,  the 
means  by  which  to  implement  their  results  is  not 
immediately  apparent. 

The  main  contributions  of  this  work  are  explicit  re¬ 
sults  for  the  probability  distribution  and  all  moments 
of  the  random  failure  time  for  a  general,  single-unit 
system  that  deteriorates  continuously  and  additively 
due  to  the  influence  of  a  random  environment  that 
is  modeled  as  a  general,  finite-state  Markov  process. 
The  model  focuses  on  the  continuous  wear  caused  by 
this  random  environment  but  ignores  any  additional 
shocks.  The  main  results  of  the  paper  are  compact, 
Laplace-Stieltjes  transform  expressions  for  the  failure 
time  distribution  and  all  moments  that  may  be  imple¬ 
mented  using  existing  numerical  inversion  algorithms. 
The  distribution  and  moments  of  the  random  failure 
time  are  extremely  important  in  reliability  analysis, 
particularly  for  the  determination  of  sound  mainte¬ 
nance  policies  in  a  variety  of  commercial,  governmen¬ 
tal,  and  military  settings. 

The  remainder  of  the  paper  is  organized  as  follows. 
The  next  section  gives  the  formal  model  description. 
In  Section  3,  we  explicitly  derive  transform  results 
for  the  failure  time  distribution  for  a  single-unit  sys¬ 
tem  subject  to  a  Markovian  random  environment.  In 
Section  4  we  derive  the  analytical  expressions  for  the 
moments  of  failure  time.  Finally,  in  Section  5,  a  nu¬ 
merical  example  is  presented  demonstrating  the  main 
results  of  the  paper. 


2.  Formal  model  description 

Consider  a  single-unit  system  subject  to  continu¬ 
ous  and  additive  deterioration  in  time  due  to  an  ex¬ 
plicit  dependence  on  the  state  of  an  external  random 
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environment.  Under  normal  operating  conditions,  the 
system  accumulates  wear  until  the  magnitude  of  its 
cumulative  wear  exceeds  a  fixed  threshold  value  x,  at 
which  time  the  system  fails.  In  particular,  the  rate  of 
deterioration  (or  wear  rate)  of  the  system  at  time  t  >  0 
is  governed  by  a  random  environment  that  is  modeled 
as  a  finite-state  Markov  process  {Z(t)\  t  ^  0}.  The 
evolution  of  the  wear  process  can  be  described  by  a 
continuous-time  stochastic  process  {X(t):  t  ^  0}  that 
assumes  values  on  the  nonnegative  real  line  which 
shall  be  denoted  throughout  by  R+.  Moreover,  we  as¬ 
sume  that  {Z(t)\  t  ^  0}  and  {X(t):  t  ^  0}  are  both 
right  continuous  and  have  left  hand  limits  everywhere. 
At  times,  for  the  sake  of  brevity,  the  processes  shall 
be  respectively  abbreviated  as  Z  and  X  with  the  un¬ 
derstanding  that  both  are  defined  on  the  parameter 
set  R+. 

The  following  definitions  are  needed  to  describe  the 
mathematical  model.  Let  Z(t)  denote  the  state  of  the 
random  environment  process  at  time  t  £  R+  and  define 
its  state  space  by  the  set  S  C  N  where  N  denotes  the 
set  of  natural  numbers.  For  this  model,  we  specifically 
assume  that  Z  has  a  finite-state  space  S  =  {1,...,AI}, 
K  £  N.  Let  R(  t )  be  defined  as  the  wear  rate  of  the  sys¬ 
tem  at  time  t  £  R+  and  define  a  nonnegative  function 
r:  S  — >  R+.  The  properties  of  the  function  r(-)  will  be 
dictated  by  the  type  of  system  under  consideration  and 
its  surrounding  environment.  Since  the  wear  rate  of 
the  system  is  explicitly  dependent  on  the  environment 
process,  the  wear  rate  process  { R(t):  t  ^  0}  assumes 
values  in  the  space  =  {r(  1 ), . . . , r(K)}. 

The  evolution  of  the  system  can  be  described  as  fol¬ 
lows.  If Z(t)=i  £  S,  then  (i)f?(t)  :=  r(Z(t))=r(i)  £  3) 
and  (ii)  the  environment  transitions  from  state  i£S 
to  state  j  £  S  at  time  t  +  e  according  to  a  Markov 
transition  function  P(/)  :=  [ Pij(t)\  where  Pij(t)  := 
P{Z(t  +  e)  =  j\Z(t)  =  /}.  The  environment  process  Z 
is  a  temporally  homogeneous,  finite-state  Markov  pro¬ 
cess  so  that  pLJ(t)  does  not  depend  on  t  for  all  i,j  £  S. 
Denote  by  X(t),  the  cumulative  wear  of  the  single-unit 
system  up  to  time  t  £  R+.  X  is  a  continuous,  additive 
functional  of  Z,  and  thus,  ( X,Z )  constitutes  a  special 
case  of  a  Markov  additive  process.  In  order  to  com¬ 
pute  the  probability  distribution  and  moments  of  the 
random  lifetime  of  the  system,  we  analyze  the  bivari¬ 
ate  process  {{X{t),Z{t))\  t  ^  0}. 

With  regard  to  the  process  {X(t):  t  ^  0},  the 
following  conditions  will  normally  be  satisfied  in 


practice: 

(i)  For  all  t,  e  5=  0,X(0)=0and X(t+e)-X(t)  ^  0, 
w.p.  1. 

(ii)  The  process  { X(t ):  f  Js  0}  satisfies  the  Markov 
property, 

P{X(t  +  e)  ^  x\ X(e)  =  y,X(u),  0  <  u  <  y} 

=  P{X(t  +  e)  ^x\X(s)  =  y}.  (1) 

In  effect,  Assumption  (i)  asserts  that  the  system  op¬ 
erates  without  intervention.  In  particular,  no  repairs 
take  place  that  can  restore  the  unit  to  a  state  of  lesser 
deterioration.  Assumption  (ii)  asserts  that  the  future 
condition  of  the  system  can  be  predicted  via  the  ac¬ 
cumulated  wear  up  to  time  t,  irrespective  of  the  his¬ 
tory  prior  to  t.  If  the  single-unit  system  has  a  known 
minimum  acceptable  wear  level  x,  it  will  fail  when 
the  cumulative  wear  exceeds  this  level.  In  this  sense, 
the  time  until  failure  for  the  unit  can  be  considered  as 
a  first  passage  time  for  the  cumulative  wear  process 
{ X(t)\  l  St  0\.  We  obtain  explicit  transform  expres¬ 
sions  for  the  failure  time  distribution,  as  well  as  all 
moments  of  this  failure  time.  In  the  following  section, 
we  first  derive  the  failure  time  probability  distribution. 

3.  Failure  time  distribution 

The  environment  process  { Z(t ):  t  ^  0}  is  assumed 
to  be  a  temporally  homogeneous  Markov  process  on 
a  finite-state  space  S  :=  { I , . . . ,  A"  } .  The  cumulative 
wear  of  the  single-unit  system  up  to  time  t  £  R+  is 
defined  by 

X(t)=  [  r(Z(u))du.  (2) 

Jo 

The  system  fails  as  soon  as  the  magnitude  of  its  accu¬ 
mulated  wear  exceeds  a  fixed  threshold  value  x.  Let  Q 
denote  the  appropriate  sample  space  for  the  Markov 
additive  process  {(X(t),Z(t)):  t  ^  0}.  For  each  sam¬ 
ple  path,  (o£Q ,  the  lifetime  of  the  system  is  given  by 
the  random  variable 

Tx(co)  =  inf^ {t:  X(t\  co)  >  x}.  (3) 

Define  the  following  distribution 

Vu(x,t)  =  P{X(t)  ^  x,Z(t )  =y|Z(0)  =  i},  (4) 
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where  V,j(x,  t )  is  the  joint  probability  that,  at  time 
t,  the  degradation  of  the  system  has  not  exceeded  a 
value  x  and  the  enviromnent  process  is  in  state  j  £  S 
given  that  the  environment  was  initially  in  state  i  £  S. 
The  objective  is  to  find  the  joint  distribution  of  ( X,Z ), 
conditional  on  i,  and  apply  the  dual  relationship  of  (3) 
to  derive  an  explicit  matrix  equation  for  the  failure 
time  distribution.  The  distribution  matrix  of  X(t)  is 
defined  as 

\{x,t)  =  [VlJ(x,t)\.  (5) 

Define  the  matrix  transforms 

\*(x,s)=  f  e~sl\{x,t)dt,  (6) 

J  R+ 

the  Laplace  transform  of  V(x,  t)  with  respect  to  t  and 

\*{u,s)  =  [  e-“V*(d x,s),  (7) 

J  R+ 

the  Laplace-Stieltjes  transform  of  \'(x,s)  with  re¬ 
spect  to  x.  Let  dij(x)  :=  VLJ(x,  0 )  and  define  a(x)  = 
[a,  y(x)  |.  It  is  clear  that 


Proof.  Fix  £  >  0. 

Vij(x,  t  +  e) 

=  P{X(t  +  £)  ^  x,Z(t  +  e)  =  jjZ(O)  =  /} 

=  ^P{Z(t  +  £)=y|X(t  +  £)^x, 

k 

Z(t)  =  k,Z(0)  =  i} 

xP{X(t  +  e)  ^  x\Z(t)  =  k,Z( 0)  =  /'} 

xP{Z(t)  =  k} 

=  (1  +  eqjjWuix  -  sr(J),t) 

+  ^  E(lkj  V'AX  ~  Er(k\ 0  +  0(e).  (12) 

kes\{j} 

Simplifying,  dividing  by  the  time  increment  s  and  let¬ 
ting  e  j  0,  it  is  seen  that  Fi j(x,t)  satisfies  the  partial 
differential  equation 

S  Vij(x,  t )  S  Vij(x,t) 

- - H - i - r(j) 


ai,Ax)  =  J{i=j}’  (8) 

where  denotes  the  indicator  function.  Further 
define  the  Laplace-Stieltjes  transform  of  a,,y(x)  with 
respect  to  x  by 

Sij(u)=  [  e-"xV,-j( du,0)  (9) 

J  R+ 

and  its  corresponding  vector  transform 

a  (u)=  j  e_'“V(dM,0).  (10) 

J  R+ 

The  double  transform  V*(w,s)  is  next  derived  un¬ 
der  the  assumption  that  the  environment  process, 
{ Z(t):  t  St  0},  is  a  finite-state  Markov  process. 


=  £^(*>0,  j  e  S  (13) 

kes 


which  may  be  written  in  matrix  form  as 


8V(x,  t) 
dt 


0V(x,O 

0X 


Rd  =  V(x,0Q- 


(14) 


Next,  we  employ  the  transform  matrix  expressions  of 
(6)  and  (7).  Taking  the  Laplace  transform  of  both 
sides  of  (14)  with  respect  to  t,  yields 

sV*(x,s)  -  a(0)  + - ;  ’  ^  Rd  =  V*(x,s)Q. 

0X 


Now,  taking  the  Laplace-Stieltjes  transform  of  the 
above  equation  with  respect  to  x  gives 


Theorem  3.1.  If  the  Z  process,  defined  on  the 
finite-state  space  S,  is  a  Markov  process  with  in¬ 
finitesimal  generator  matrix,  Q  =  \q, j\,  i,j  £  S,  then 


s\*(u,s )  -  I  +  mV*(m,,s')R«  =  V*(u,s) Q.  (15) 

The  theorem  follows  by  rearranging  the  terms  of  ( 1 5 ). 

□ 


V*(m,5)  =  [uRZ3+5I-Q]^1, 

Re(u)>  0,  Re(s)  >  0.  (11) 

where  R d  =  diag(r(  \ ), . . .  ,r(K)). 


Theorem  3.1  may  now  be  used  to  derive  the 
explicit  expression  for  failure  time  distribution  of  the 
single-unit  system  subject  to  a  Markovian  environ¬ 
ment.  Define  the  joint  distribution  of  Tx  and  Z(t), 
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conditional  on  Z( 0),  by 

W,Jx,t)  =  P{Tx  ^  t,Z(t)  =j\Z(0)  =  i}. 

Due  to  the  dual  relationship  of  (3),  we  obtain  the 
marginal  distribution  of  Tx,  conditioned  upon  i  £  S  as 

W,(x,t)=  l-22vu(x>0-  (16) 

jes 

It  follows  that  G(x,  t )  :=  P{TX  ^  t},  the  unconditional 
distribution  of  Tx,  is  given  by 


conditioned  upon  the  event  that  the  initial  state  of  the 
environment  Z  is  i  £  S.  The  following  lemma  will  be 
needed  to  derive  an  explicit  expression  for  the  (condi¬ 
tional)  moments  of  the  random  time  to  failure  of  the 
system. 

Lemma  4.1.  For  n  Js  0,  the  nth  order  partial  deriva¬ 
tive  of\*(u,s )  ii'.  r.  t.  s  is  given  by 

V\u,s)  =  (-iyn\(uRD  +sl  -  Qr”"1. 
osn 


G(x,  t)  =  1  -  «V(x,  f)l,  (17) 

where  a  :=  [a,]  with  a,  :=  P{Z(0)  =  /}  and  1  is  a 
K -dimensional  column  vector  of  ones.  We  now  define 
the  following  transforms, 

G*(x,s)=  I  e~s‘G(x,t)dt, 

G*(u,s)  =  [  e~wcG*(dx,s). 

J  R+ 

Theorem  3.2.  The  failure  time  distribution  of  a 
single-unit  system  in  a  Markovian  environment  is 

G*(u,s)= - aV*(M,s,)L  Re(u)>  0,  Re{s)>  0. 

(18) 


Proof.  The  proof  follows  directly  by  taking  the 
Laplace  transform  of  (17)  with  respect  to  t  and  then 
taking  the  Laplace-Stieltjes  transform  with  respect 
to  x.  □ 

It  should  be  noted  here  that  the  two-dimensional 
transform  of  (18)  can  be  inverted  numerically  using 
algorithms  such  as  those  due  to  Abate  et  al.  [1]  and 
Moorthy  [9].  In  Section  4,  we  use  the  results  of  The¬ 
orem  3.1  to  derive  transfonn  expression  for  all  mo¬ 
ments  of  the  random  lifetime  Tx. 


4.  Analysis  of  moments 

We  now  turn  our  attention  to  the  moments  of  Tx 
whenever  the  threshold  value  x  £  R+  is  finite.  Define 
the  conditional  expectation  of  the  nth  moment  of  Tx  by 

f!(x)  =  E[Txn\Z(0)  =  i]  (19) 


Proof.  The  lemma  will  be  proved  by  induction  on  n. 
By  (11),  the  result  clearly  holds  for  n  =  0.  Let  n  =  1 
and  note  that 


V*(m,5)(mRd  +  si  —  Q)  =  I 

also  by  (11).  Differentiating  both  sides  of  the  above 
equation  with  respect  to  s  shows  that 

—Y*(u,s)  =  (-1)Y(u,s)(uRd  +sl  —  Q)_1, 

05 

so  the  result  holds  for  n  =  1.  For  the  inductive  step, 
assume  that 

0m 

—  Y*(u,s)  =  (-l)mm\(uRD  +sl  -  Q)-'"-1  (20) 

csm 

for  an  arbitrary  m  £  N.  Rewriting  (20)  gives, 

- 

—  Y*(u,s)(uRD  +  si  -  Q)m+1  =  (-ir«!  x  I. 
osm 

Differentiating  both  sides  of  the  above  equation  and 
rearranging  terms,  we  have 

0m+l 

—  V>,5)(MRZ3+.vI-Qr+1 


=  (-l  )(m  +  1)^  V>,5)(uRd  +  si  -  Q)m 

Substituting  (20)  into  the  above  equation,  we  obtain 

0'H+1  ~ 

— - — r  Y*(u,s) 

dsm+ 1  v  ’ 


=(-l  )'n+\m  +  1)!(«R0  +sl  -  Q)-"‘-2 
and  the  proof  is  complete.  □ 


We  next  state  our  result  for  the  conditional  nth  mo¬ 
ment  (n  ^  1 )  of  the  failure  time. 
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Theorem  4.1.  For  n  >  1,  the  Laplace— Stieltjes 
transform  of  £"(x)  is 

|"(w)  =  n\(uRD  —  Q)_11.  (21) 

Proof.  By  Eq.  (16),  the  Laplace-Stieltjes  transform 
of  Wj(x,  t)  with  respect  to  t  may  be  written  as 

Wi{x,s)=\-Y,sK)M,  (22) 

j&S 

where  Vf(x,s)  is  the  Laplace  transform  of  Vjj(x,t ) 
with  respect  to  t.  It  is  well-known  that  the  conditional 
expectation  of  (19)  may  be  obtained  by  evaluating  the 
nth  order  derivative  of  (22)  at  s  =  0.  To  that  end,  let 

' lf(x,s) 

:=(-iy^-lwi(x,s)  =  (-iy+l 

V-  0"  *  0"_1  *  \ 

I]5  g Jn  VuM>s)  +  n  Vu(x, s)  .  (23 ) 

jes  / 

Eq.  (23 )  implies  that 


applying  Lemma  4.1  and  setting  s  =  0,  we  have  that 

£(«)  =  nKuRfl-Qr'l 

and  the  proof  is  complete.  □ 

Let  c"(x)  :=  E\T"]  be  the  unconditional  nth 
moment  of  the  failure  time.  The  following  corollary 
follows  directly  from  Theorem  4.1. 

Corollary  4.1.  Suppose  Z  has  initial  probability  dis¬ 
tribution  a.  Then, 

<f(w)  =  n!a(MRD-QrIl.  (26) 

The  above  transform  can  be  easily  inverted  numeri¬ 
cally  by  a  number  of  one-dimensional  inversion  algo¬ 
rithms.  The  reader  is  referred  to  Abate  and  Whitt  [2] 
for  one  such  approach.  In  the  following  section,  the 
main  results  of  the  paper  will  be  demonstrated  on  an 
example  problem. 


5.  Numerical  example 


%(X)=r(x,  o) 


=(->r‘»  Ei 


jes 


(24) 


s= 0. 


In  order  to  solve  the  differential  equation  (24),  trans¬ 
form  methods  are  again  employed.  The  Laplace- 
Stieltjes  transform  of  f(x)  with  respect  to  x  is 

ti(u)=  f  e-“S(dx). 

J  R+ 


By  taking  the  Laplace-Stieltjes  Transform  of  Eq.  (24) 
on  both  sides, 


an- 1 


jes 

0„-i 


i=0 


=  (^1} 


1, 


(25) 


s= 0 


where  \*(u,s )  is  the  matrix  transform  of  (11)  and  1 
is  a  -dimensional  column  vector  of  ones.  Finally, 


In  this  section,  we  demonstrate  the  applicability  of 
the  main  results  of  this  paper  via  a  numerical  example. 
In  particular,  we  consider  the  propagation  of  a  crack 
in  metallic  materials.  For  a  more  thorough  treatment 
of  fatigue  crack  dynamics,  the  reader  is  referred  to  the 
paper  by  Ray  and  Tangirala  [12]  and  references  con¬ 
tained  therein.  Assume  that  the  metallic  component 
is  placed  in  operation  in  perfect  working  order  (i.e., 
there  is  no  crack  in  the  material  initially).  However, 
due  to  normal  wear  and  fatigue,  a  crack  is  initiated 
and  continues  to  grow  in  time. 

Let  X(t)  denoted  the  length  of  the  crack  at  time  t  and 
assume  that  the  (linear)  rate  at  which  the  crack  grows 
is  subject  to  its  random  environment  (applied  stress, 
ambient  conditions,  and  other  factors).  It  is  assumed 
that  these  environmental  factors  can  be  characterized 
by  {Z(0:  t  Js  0},  a  temporally  homogeneous  Markov 
process  that  alternates  between  two  distinct  states  such 
that  its  finite-state  space  is  given  by  S  =  {1,2}  (i.e., 
Z  is  a  simple  on/off  process).  Whenever  the  environ¬ 
ment  is  in  state  i  £  S,  the  crack  grows  at  rate  r(i)  units 
per  unit  time,  i  £  S.  The  initial  distribution  of  the  envi¬ 
ronment  process  is  arbitrarily  chosen  to  be  a  =  [1  0], 
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Table  1 

P{TX  ^  l }  when  Z  is  an  on/off  process 


t 

Analytical 

Simulation 

t 

Analytical 

Simulation 

1.197 

1.2589E-01 

1.2880E-01 

2.018 

9.6765E-01 

9.6980E-01 

1.288 

2.3732E-01 

2.4300E-01 

2.109 

9.8299E-01 

9.8290E-01 

1.379 

3.7196E-01 

3.7040E-01 

2.200 

9.9149E-01 

9.9140E-01 

1.470 

5.1284E-01 

5.1400E-01 

2.291 

9.9584E-01 

9.9570E-01 

1.562 

6.4366E-01 

6.4840E-01 

2.382 

9.9819E-01 

9.9800E-01 

1.653 

7.5388E-01 

7.5010E-01 

2.474 

9.9912E-01 

9.9940E-01 

1.744 

8.3960E-01 

8.3600E-01 

2.565 

9.9950E-01 

9.9960E-01 

1.835 

9.0102E-01 

9.0020E-01 

2.656 

9.9989E-01 

9.9960E-01 

1.926 

9.4205E-01 

9.4240E-01 

2.747 

9.9997E-01 

1.0000E+00 

i.e.,  the  environment  starts  in  state  1  with  probabil¬ 
ity  1 .  Moreover,  the  infinitesimal  generator  matrix  is 
given  by 


Q  = 


-b 


b 


while  the  diagonal  matrix  of  wear  rates  is 
r(l)  0 


R 


0  r(  2) 


Now,  the  partial  differential  equation  describing  the 
probability  distribution  of X(t)  is  given  by 

£  Vij{x,  t)  dVij(x,t) 


dt 


0X 


r(j) 


2 

=  qkJ  V‘AX’  0,  7  =  1,2 

k=  1 


which  can  be  written  in  matrix  form  as 


8V(x,0 

dt 


d\(x,t) 

0X 


R/)  =  V(x,r)Q, 


where  V(x,t)  =  \VLJ(x,t)\.  Now  applying  Eqs.  (11) 
and  (18),  the  transform  G*(u,s )  is  obtained  by 


G*(u,s)  =  -  —  a(nR0  +  si  —  Q)  'l 
s 

1 
1 


=  -  —  [1  0](<5y  —  ab) 


-l 


<5  b 

a  y 


where  <5  =  r(2)u  +  s  +  a  and  y  =  r(l)u  +  s  +  b 
are  scalars.  The  above  two-dimensional  LST  is  con¬ 
verted  to  a  two-dimensional  Laplace  transform  by 


Table  2 

Comparison  of  £*(x)  when  Z  is  an  on/off  process 


X 

Analytical 

Simulated 

0.25 

3.6094E-01 

3.6139E-01 

0.50 

7.3594E-01 

7.3617E-01 

1.00 

1.4859E+00 

1.4852E+00 

5.00 

7.4859E+00 

7.4875E+00 

10.00 

1.4986E+01 

1.4979E+01 

20.00 

2.9986E+01 

2.9987E+01 

50.00 

7.4986E+01 

7.4983E+01 

pre-multiplying  both  sides  of  the  equation  by  u~l .  The 
inversion  algorithm  of  Moorthy  [9]  was  implemented 
to  compute  the  analytical  cumulative  probability  val¬ 
ues  at  several  points  and  compared  with  simulation 
results. 

The  specific  problem  parameters  chosen  for  this 
example  are  as  follows.  The  threshold  value  x  =  1.0 
and  the  values  comprising  the  generator  matrix  are 
a  =  b  =  25/3.  The  state-dependent  crack  growth  rates 
are  r(l)=  1.0833  and  r(2)=0.250.  For  comparison 
purposes,  an  empirical  cumulative  distribution  func¬ 
tion  was  constructed  via  Monte  Carlo  simulation  tech¬ 
niques.  The  cumulative  probabilities  were  compared 
at  several  values  of  t  and  the  results  summarized  in 
Table  1. 

In  the  following  example,  we  demonstrate  com¬ 
putation  of  the  first  two  moments  of  the  failure  time 
using  identical  problem  parameters.  In  this  case,  the 
one-dimensional  inversion  algorithm  of  Abate  and 
Whitt  [2]  was  implemented  and  found  to  perform 
extremely  well.  Tables  2  and  3  contain  a  summary  of 


244 


J.P.  Kharouf eh  I  Operations  Research  Letters  31  (2003)  237-244 


Table  3 

Comparison  of  i2(x)  when  Z  is  an  on/off  process 


X 

Analytical 

Simulated 

0.25 

1.4462E-01 

1.4495E-01 

0.50 

5.7353E-01 

5.7360E-01 

1.00 

2.275 1E+00 

2.275 1E+00 

5.00 

5.6388E+01 

5.6407E+01 

10.00 

2.2528E+02 

2.2537E+02 

20.00 

9.0056E+02 

9.0038E+02 

50.00 

5.6264E+03 

5.6260E+03 

the  numerical  values  as  compared  with  Monte  Carlo 
simulation  results. 
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